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THEORETICAL  PRESSURE  DISTRIBUTIONS  FOR  SOME  SLENDER 
WING-BODY  COMBINATIONS  AT  ZERO  LIFT 
By  Paul  F.  Byrd 


SUMMARY 


Pressure  distributions  are  calculated  for  some  symmetrical  wing-body 
combinations  at  zero  lift.  The  theory  of  the  calculations  is  based  on 
the  assumption  of  extremely  slender  wings  and  bodies  and  yields  results 
for  both  subsonic  and  supersonic  speeds.  The  examples  considered  are 
swept  wings  of  constant  chord  mounted  on  bodies  of  nearly  cylindrical 
form. 


Of  particular  interest  is  the  effect  of  indenting  the  body  on  the 
distribution  of  pressure  over  the  wing.  When  the  indentation  is  such  as 
to  maintain  a  constant  total  area  of  the  cross  sections  normal  to  the 
stream,  the  theoretical  pressure  disturbances  remain  small  throughout  the 
transonic  range.  With  such  indentation  the  isobars  tend  to  remain  smooth 
and  nearly  parallel  to  the  sweep  of  the  wing  surface. 


INTRODUCTION 


In  several  papers,  important  extensions  to  the  Munk-Jones  slender- 
body  theory  (refs.  1  and  2)  for  lifting  wings  and  bodies  have  been  made 
to  include  the  theoretical  effects  of  thickness  on  the  aerodynamics  of 
wings  and  wing-body  combinations.  Ward  (ref.  3) >  solving  the  linearized 
differential  equation  for  the  perturbation  velocity  potential  by  opera¬ 
tional  methods,  and  employing  asymptotic  expansion  of  the  solution, 
investigated  the  flow  around  bodies  of  general  cross  section  at  supersonic 
flight  speeds.  By  a  different  procedure,  similar  results  for  a  wing, 
body,  or  wing-body  combination  at  subsonic  speeds  have  been  developed 
by  Heaslet  and  Lomax  (refs.  4  and  5).  An  analysis  for  subsonic  flow  was 
also  carried  out  independently  by  Adams  and  Sears  (ref.  6)  who,  in  addi¬ 
tion,  made  an  extension  for  not -so -slender  wings.  Confining  themselves 
to  wing6  at  zero  angle  of  attack,  Keune  (ref.  7)  and  Oswatitsch  and  Keune 
(ref.  8)  have  recently  obtained  a  slender-body  theory  that  is  slightly 
different  from  those  of  references  4,  and  6.  In  a  more  recent  report. 
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Harder  and  Klunker  (ref.  9)  have  applied  the  basic  ideas  of  the  slender- 
body  approximation  to  the  nonlinear  transonic  equation  for  the  velocity 
potential. 

The  principal  object  of  the  present  investigation  is  to  apply  the 
general  method  of  reference  1+  in  calculating  the  pressure  distribution 
for  some  special  cases  of  nonlifting  slender  wing-body  combinations  in 
subsonic  and  supersonic  flow.  The  wing  of  the  combinations  is  swept  back 
and  has  a  symmetrical  section  with  rounded  leading  edges.  Determination 
is  made  of  the  pressure  for  the  wing  alone  and  for  cases  when  the  wing  is 
mounted  on  a  circular  cylinder  or  combined  with  a  body  indented  such  that 
the  axial  variation  of  cross-sectional  area  of  the  combination  is  constant. 
The  effects  of  Mach  number  and  sweep  angle  are  included  in  the  results 
presented. 


LIST  OF  IMPORTANT  SYMBOLS 


C*r, 


If 


m 

M 

r 

Re 

Ro 

R(x) 

sgn(x-f) 

so 


value  of  x  at  which  sQ  =  s(x) 
root  chord 


pressure  coefficient, 


-2u 

Un 


pressure  coefficient  on  indented  wing-body  combination 
over-all  length  of  the  wing 
slope  of  wing  leading  edge  (See  sketch  (f).) 
free -stream  Mach  number 
polar  distance  in  y,z  plane 
real  part  of  a  complex  quantity 
radius  of  cylindrical  portion  of  body 
radius  of  indented  body  of  revolution 
sign  of  (x-$) 
maximum  value  of  s(x) 


(/y2  +  z2^ 


s(x) 

S(x) 


local  semispan 

local  cross -sectional  area  of  wing  alone 
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t(x)  spanwise  distance  from  x  axis  to  wing  trailing  edge 

u  perturbation  velocity  in  the  x  direction 

Uo  velocity  of  free  stream 

v  perturbation  velocity  in  the  y  direction 

vr  radial  component  of  velocity  in  yz  plane 

w  perturbation  velocity  in  the  z  direction 

x,y, z  Cartesian  coordinates  (x  downstream,  y  to  starboard,  z 

upward) 


Q  polar  angle  in  yz  plane 

X(x,y)  slope  of  wing  surface  in  x  direction 

<p  perturbation  velocity  potential 

cp2  part  of  potential  satisfying  cpyy  +  <pzz  =  0 

£  complex  variable  (y  +  iz) 

complex  variable  (yi  +  izi) 
t  constant  related  to  tq 

tq  maximum  thickness  of  wing  section 

Subscripts 

B  body 

l  lower  surface  of  the  wing  (z  =  0  plane) 

sub  subsonic 

8 up  supersonic 

u  upper  surface  of  the  wing  (z  ■  0  plane) 


W 


wing 
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ANALYSIS 


The  thickness  distribution  of  a  symmetrical  wing  is  prescribed. 

Under  the  restriction  that  the  thickness  is  small  and  that  the  configura¬ 
tion  is  slender,  formulas  to  be  applied  later  (see  "Applications"  section) 
will  now  be  briefly  presented  for  determining  the  pressure  coefficients 
for  particular  nonlifting  wing-body  combinations.  Equations  are  first 
given  for  the  wing  without  a  body,  and  are  then  modified  for  cases  when 
the  wing  is  mounted  on  a  circular  cylinder  or  on  an  indented  body  of 
revolution. 


Nonlifting  Wing  Alone 


The  differential  equation  and  boundary  conditions.-  Expressed  in 
terms  of  the  perturbation  potential  cp(x,y,  z),  the  basic  linearized  par¬ 
tial  differential  equation  for  subsonic  as  well  as  supersonic  flow  is  the 
familiar  Prandtl-Glauert  equation 


(i  -  M02)9XX  +  <Pyy  +  «P22  =  0  (1) 

where  M0  is  the  free-stream  Mach  number.  If  the  surface  of  a  wing 
z(x,y)  is  given,  solutions  to  the  differential  equation  must  satisfy  the 
boundary  condition  that  the  flow  is  parallel  to  the  wing  surface.  When 
the  wing  is  thin,  it  is  sufficient  to  satisfy  this  requirement  in  the 
plane  z  =  0.  Analytically,  the  expression  of  the  condition  is 


wu(x,y) 


Uo 


dz 


u 


3x 


UoMx>y) 


(2) 


where  U0  is  the  free-stream  velocity  and  wu(x,y)  is  the  vertical 
induced  velocity  on  the  upper  side  of  the  z  *  0  plane. 

Velocity  potential.-  When  the  flow  is  supersonic,  the  formula  for 
the  perturbation  potential  subject  to  the  boundary  condition  (2)  is  known 
to  be  (e.g.,  see  ref.  10) 


- an  (3) 

3 *J (y  -  n)2  +  z2 


cp(x*y>z) 


arccosh 
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where  the  region  of  integration  T  is  the  portion  of  the  plan  form  lying 
within  the  Mach  forecone  from  the  point  x,y, z.  In  the  case  of  subsonic 
flow  fields,  the  solution  of  the  differential  equation  (l)  may  be  expressed 
in  the  form 


<P(x,y,z) 


vu(M) 


arcsinh 


— ■  — - d£  dq 

p/(y  "  T))2  +  z2 


(b) 


with  T  now  extending  over  the  entire  plan  form.  These  two  solutions 
yield  the  potential  due  to  a  distribution  of  sources  of  strength  propor¬ 
tional  to  the  slope  of  the  wing  surface  \j(x,y). 


When  the  wing  is  slender,  that  is,  if 


£ 


A 


2  +  z2 


is  considered 


x  -  i 

very  small,  further  approximations  to  the  linearized  potential  (3)  and  (4) 
can  be  readily  made.  In  this  event,  one  may  employ  the  approximate 
relations 


arccosh 


x  -  l 

3^/(y  -  tj)2  +  z2 


~  In 


2(x  -  |) 

3-/(y  -  t)2  +  z2 


>  (5) 


arcsinh 


x  -  § 

P^Ay  -  n)2  + 


=  *  sgn(x  -  5)ln 
z2 


2|x  -  M 

3</(y  -  n)2  +  z2 


j 


where  the  symbol 
taken. 


sgn(x  -  5)  means  that  the  sign  of  (x  -  I)  is  to  be 


Consider  now  a  thin  pointed  wing 
of  symmetrical  section  with  straight 
or  sweptforward  trailing  edges  as  is 
shown  in  sketch  (a).  Use  of  the 
relations  (5)  in  equations  (3)  and  (4) 
then  gives  the  result  (ref.  4) 


<p(x,y,z)  =  <p2(x,y,z)  +  g(x) 


(6) 
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where 

U  ^a(x) 

92(x,y,z)  =  — J  Xu(x,T))ln  [(y  -  n)2  +  z2]dn  (7) 

s2(x) 

and 

g{x)='^hf  S'(6)ln  d&  (8a) 


for  supersonic  flow,  and 

Uo  8  r°  2  lx  -  i  I 

g(x)  *  -  ujj-  ^ J  sgn(x  -  l)S»(|)ln  - - - -  d|  (8b) 


for  subsonic  flow.  The  function  S’(x),  the  derivative  of  the  cross- 
sectional  area  of  the  wing  in  a  yz  plane,  is  found  from 

sa(x) 

S»(x)  =  2  zu(x,n)dT)  (9) 

s2(x) 

It  is  seen  from  equation  (6)  that  the  slender-body  approximation  to 
the  linearized  supersonic  and  subsonic  potentials  consists  of  two  parts. 
The  first  part  ($2)  in  each  case  is  independent  of  Mach  number  and  is  a 
harmonic  function  in  the  transverse  plane;  that  is,  it  satisfies  the  two- 
dimensional  Laplace  equation 


‘Pyy  +  ‘Pzz  =  0 


(10) 


The  second  part  (g)  depends  on  the  cross-sectional  area  of  the  wing  and 
is  a  function  of  x  and  3  only# 

Inspection  of  equations  (6)  and  (8)  shows  that  the  value  of  the 
potential  for  a  particular  wing  at  some  Mach  number  Mi  can  be  written 
in  terms  of  the  potential  given  at  another  Mach  number  M0.  Thus 
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Ks^P3l 


^  ♦  sS  s.Wi»  a 


(lla) 


Uo  Si 

1>8ubeo+  2^  S-(x)ln  — 


with  pi  =  y|Mi2  -  l|  and  p0  =  J  |M02  -  l| .  The  supersonic  and  subsonic 
potentials  may  also  be  related.  If  one  has  already  found  an  expression 
for  the  supersonic  potential,  he  can  then  obtain  the  subsonic  potential 
by  means  of  the  equation 


^sub 


+  *7  S»(x)ln 


u0  r 
W 


s«  (Q 

x  -  i 


(lib) 


-y=s(x) 


Certain  symmetrical  wings  whose  trailing  edges  are  swept  back  as  in 
sketch  (b)  can  also  be  treated  by  the  simplified  theory,  provided  the 
chordwise  variation  in  shape  of  the 
cross  sections  is  sufficiently  smooth  J Uo 
and  gradual  that  the  assumption  of 

two-dimensional  transverse  flow  may  q _ „ 

reasonably  be  applied.  Since  the  n  t\ 

wing  is  at  zero  angle  of  attack  and  /  \ 

the  flow  is  symmetrical  about  the  " y  V'  ‘  ' 

xy  plane,  the  source  distribution  0 

Au(x,y)  in  the  plane  of  the  wing  in  /  '  \  U 

region  W  between  the  axis  and  the  /  A  ° 

trailing  edge  is  set  equal  to  zero.  y=-t(x)  /  ~P\  '-y=t(x) 

If  the  limits  of  integration  are  V/ 

properly  adjusted,  equation  (6)  will 

then  formally  still  apply  and  yield  Is  ^ 

expressions  for  the  potential  in  the  ’  1 

two  regions  1  and  2  of  the  sketch.  * 

Thus  for  region  1,  in  the  z  *  0  plane, 

one  obtains  Sketch  (b) 


Sketch  (b) 


'psup  =  %  f  xu(x,n)in|y*  -  T»eldn  • 


x 

S  hi  Sl,(6)ln  f  (x  "  *)d*'  (x^co)  U2a) 
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and  in  region  2 


Uc  rs  U0  £  So 

‘Psup  =  T  J  ln  |y2  ■  T)2  |dTl  “  2 it  si'(£)ln(x  -  l)d6  - 

r\ 


S2'(0ln  |  (x  -  (x  >  C0) 


(12b) 


where 


s(x) 


si'(x)  =  4  ~ J*  zu(x,n)dT),  0  <  X  <  cQ 


s(x) 


Sa'U)  4  Sc  /  ^u(x,n)dn,  cQ  <  x  <  ZG  J 
t(x) 


(13) 


The  subsonic  solution  is  obtained  from  the  supersonic  solution  by  using 
equation  (lib). 

Pressure  coefficient.-  After  determining  the  potential  for  a  wing 
from  the  equations  in  the  foregoing  section,  the  pressure  on  the  surface 
of  the  wing  is  found  by  differentiation.  The  pressure  coefficient  is 
related  to  the  perturbation  velocities  by  the  equation  (ref.  4) 


cp(x,y,0) 


2d©  1 

L5oa* +  57 


1  f  dq>Y 

07  WJ 


z=0 


(1W) 


which  is  an  approximation  to  the  complete  Bernoulli  equation  consistent 
with  both  the  linearized  differential  equation  (1)  and  the  assumption  of 
slenderness.  For  planar  problems,  further  simplification  achieved  by 
neglecting  the  nonlinear  terms  (dq)/dy)2  and  (dqp/dz)2  yields  satisfactory 
estimates  for  the  coefficient. 
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Thus 


cp(x,y,0) 


(l4b) 


where  9  is  obtained  from  equations  (lib)  and  (12). 


It  should  be  pointed  out,  however,  that  equations  (lib),  (12),  and 
(l4)  will  not  furnish  realistic  results  for  the  pressures1  for  a  wing  in 
either  subsonic  or  supersonic  flow  unless  certain  restrictions  are  imposed 
on  the  gradient  of  cross-sectional  area  and  its  derivative  S"(x).  Con¬ 
sider,  for  example,  the  pressure  found  along  the  line  AB  where  regions 
1  and  2  in  sketch  (b)  join.  In  supersonic  flow,  one  finds  from  equa¬ 
tions  (12)  and  (l4b)  that  at  this  line  there  is  a  jump  Acp  in  the  pres¬ 
sure  coefficient  given  by 


Acp  -  -  <*.  ■ ;  »7.*u0') 


4 


Au(x»t) 


In 

c0  +  e 


(15) 


where 


ASn( c0)  =  S!"(c0)  -  S2"(c0) 


and  where  the  usual  assumption  that  S'(0)  =  Sf(c0)  =  0  has  been  made. 
Evidently,  along  the  line  AB  a  logarithmic  singularity  will  occur  in 
the  pressure  distribution  in  going  from  region  1  to  region  2  if 
Si"(co)  /  S2M(c0).  (The  infinity  would  of  course  be  higher  than  loga¬ 
rithmic  if  AS"(c0)  is  singular.)  The  formulas  for  the  pressure  coef¬ 
ficients  presented  in  this  report  and  in  reference  4  are  therefore  good 
only  for  cases  where  the  plan  form  and  the  slope  of  the  wing  surface  are 
sufficiently  smooth  so  that  there  are  no  abrupt  changes  in  either  S'  or 
Sn.  (This  holds  true  also  for  the  formulas  given  in  references  3,  6,  and 
7.  The  restriction  on  S",  however,  may  be  somewhat  relaxed  in  employing 
the  slender-body  theory  of  reference  8.) _ 

Although  the  theory  may  give  spurious  infinite  pressures  on  certain 
portions  of  the  wing,  the  results  obtained  for  the  wave  drag  by  integrat¬ 
ing  the  product  of  pressure  and  surface  slope  over  the  wing  may  be  finite 
and  reasonable. _ _ 
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Even  if  there  is  no  Jump  in  the  pressure  coefficient  at  the  line  AB, 
the  slope  of  the  pressure  curve  may  not  be  continuous.  It  can  be  shown, 
for  instance,  that  a  singularity  will  occur  in  the  slope  of  the  curve  for 
the  pi’essure  coefficient  in  case  Sm  (c0)  is  discontinuous. 


Wing  on  a  Cylindrical  Body 


The  equations  presented  in  the  foregoing  sections  for  the  wing  alone 
will  now  be  modified  to  yield  formulas  for  calculating  the  pressure  dis¬ 
tribution  on  a  combination  composed  of  a 
I  (Jo  symmetrical  wing  zu(x,y)  =  -Zj(x,y) 

mounted  on  an  infinite  circular  cylinder 
having  radius  R0(see  sketch  (c)).  The 
surface  of  the  wing  chosen  here  will  also 
be  considered  symmetrical  about  the  xz 
plane,  but  the  method  applies  equally  well 
if  this  is  not  the  case.  The  procedure 
followed  is  essentially  the  same  if  the 
fuselage  is  any  body  of  revolution  instead 
of  a  circular  cylinder. 

In  studying  such  combinations,  it  is 
usually  convenient  to  introduce  a  second 
coordinate  system.  Let  the  yz  plane  be 
represented  by  a  complex  variable 

C  =  y  +  iz  =  reie  (16) 

and  then  consider  a  Ci  plane 
Ci  =  yi  +  izi  =  riei0i 

obtained  from  the  C  plane  by  the 
Joukowski  transformation 
Sketch  (c)  2 

Cx  =  C  +  (17) 

The  transformation  (17)  maps  the  C  plane  onto  the  Ci  plane  so 
that  a  circle  representing  a  section  of  the  body  in  the  £  plane  is 
mapped  onto  a  portion  of  the  real  axis  in  the  Ci  plane,  while  the  part 
of  the  real  axis  outside  the  circle  is  transformed  into  an  adjoining  part 
of  the  real  axis  of  the  £*  plane.  (See  sketch  ( c) . )  It  can  easily  be 
shown  from  the  equation  that  the  geometric  relations 
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Rx  =  2R0 

Ro2 

yi  -  y  +  -j-  , 

yi  =  2R0cos  0, 
R02 

si  =  s  +  -a- 


(18) 


hold  for  Zi  equal  to  zero. 

The  technique  employed  is  to  transform  the  boundary  conditions  from 
the  £  plane  to  the  plane  and  to  find  a  solution  in  the  latter 

plane.  The  solution  is  then  transformed  back  to  the  physical  plane  for 
the  completion  of  the  problem. 

Perturbation  velocity  potential  at  the  wing  surface.-  A  consequence 
of  the  conformal  transformation  is  that  the  complex  velocities  in  the 
two  planes  £  and  are  related  by  the  equation 


From  these  equations  it  follows  that  the  boundary  conditions  in  the  zi  =  0 
plane  are 
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w  »  vr  =  U0(dR0/dx) 
1  "  2  sin  0  2  sin  0 


yi2  <  R12 


>  (20) 


Vi 


R12  <  yi2  <  Si2 


with  w  related  to  the  streamwise  slope  of  the  wing  surface  by  equa¬ 
tion  (2) . 


The  two-dimensional  solution  in  the  £1  plane  at  zi  =  0  is, 

from  equation  (7)> 


Putting  now 


<p2(x»yi»°)  =  I J  wi(x>niUn|yi  -  ni km 

-81 

Ro2 

m  -  n  +  — 


(21) 


and  using  relations  (18)  and  the  boundary  conditions  (20),  one  finally 
has  for  equation  (21) 


tt  8(x) 

U0  r 

92(**y)  =«“/ 


(y2  -  n2)(y2n2  -  R04) 
yV 


dq 


(22) 


If  one  of  the  functions  g  given  by  formulas  (8)  is  then  added  to  this 
equation  for  <Pg ,  the  velocity  potential  on  the  wing  may  be  obtained  for 
either  supersonic  or  subsonic  flow.  Thus,  for  the  supersonic  perturbation 
potential  on  the  upper  surface,  there  results 


s(x) 

<Psup  =  icj  ^u(x>Tl) ln 
Ro 


(y2  -  n2)(y*n8  -  Ro4) 

yV 


H°±  /  S’(5)ln|  (x  -  |)d6  (23) 

2«  dxg  3 


and  for  the  subsonic  potential 
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The  coefficient  for  a  particular  wing  at  some  Mach  number  Mi,  written 
in  terms  of  the  coefficient  at  another  Mach  number  M0,  is 


♦  ;  I 


(27a) 


-P  k  =  CP  v  +  n  r 

t^subj31  ^subj30  71  pi 


(27b) 


where  =  J  jMi2  -  l|  and  fi0  =  J  |m02  -  1 1 . 

It  should  be  noted  that  formulas  (25a)  and  (25b)  are  not  uniformly 
valid  for  all  values  of  R0.  If  the  radius  R0  is  equal  to  zero,  the 

(equations  reduce  immediately  to  those 
Uo  given  previously  for  the  wing  without 

a  body.  As  the  radius  approaches 
infinity,  however,  one  finds  that  for¬ 
mulas  (25a)  and  (25b)  yield  a  value  of 
the  pressure  coefficient  that  differ 
from  the  results  for  the  wing  alone  by 


T  £$ 


a  term  equal  to 


-  In  i  •  Employing 


Sketch  (d) 


the  equations  for  values  of  R0  that 
are  large  violates  the  assumption  of 
slenderness  according  to  the  approximate 
theory  used  herein. 

When  the  wing  mounted  on  the  com¬ 
bination  is  swept  back,  as  shown  in 
sketch  (d),  the  slope  Au(x>n)  is  taken 
equal  to  zero  in  the  gap  between  the 
trailing  edge  and  the  body.  The  for¬ 
mulas  obtained  in  regions  1,  2,  and  3 
for  the  supersonic  pressure  coefficient 
on  the  upper  surface  of  the  wing  are 
then  as  follows: 


Region  1: 


2  a_  r  .  ,  (y2  -  i2)(yV  -  R04)  .  + 

'"W  Xu(x,n)ln  yV 

Ho 

7  J  Si  *  (5) in  |  (x  -  5)d5,  (0  <  x  £  c0) 
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Region  2: 

s(x) 

cf2  -  - 1 !;/ 

t(x) 


(y2  -  n2)(y2ng  -  Ro4) 

y^r)2 


dtj  + 


rP O 


f  si'(Oln(*  -  ft )dft  +  “  a 


*  8xs 


Sa'dJln  |  (x  -  ft  )d.ft 


Co 


3 


(cQ  <  x  <  bQ)  (28b) 


Region  3: 


so 


-P3  =  -  I  £  /  *u<x»H>ln 
t(x) 


(y2  -  t)2) (yV  -  R04) 

A2 


drj  + 


I  £5/  si*(6)ln(x  -  4)d£  +  |  ~s  f  S2»(0ln(x  -  Od6  + 
o  Co 


x 

7  £l  / S3»(0ln  |  (x  -  (b0  $  x  <  lQ) 


(28c) 


where  b0  is  the  value  of  x  at  which  sQ  =  s(x).  Use  of  the  relation 
(26)  in  conjunction  with  the  above  three  formulas  will  furnish  equations 
for  the  pressure  coefficient  in  subsonic  flow  in  the  various  regions. 


Wing  on  an  Indented  Body- 


In  the  previous  formulas  for  the  pressure  distribution  attributable 
to  thickness,  the  coefficient  becomes  infinite3  when  the  Mach  number  g 
approaches  unity  because  the  formulas  contain  a  term  involving  S'(x)ln  — . 

It  is  thus  possible  to  construct  slender  configurations  which  will  give 
a  theoretically  finite  pressure  coefficient  even  at  the  speed  of  sound. 
Under  the  assumptions  of  slender-body  theory,  this  may  be  achieved  if  the 
gradient  of  cross-sectional  area  of  the  configuration  in  a  yz  plane 
vanishes  identically.  Combinations  constructed  in  this  manner  on  an 

- g - — . — - - - - -  - -  - —  ~ 

This  of  course  is  not  a  property  of  slender-body  theory  alone. 

Except  in  particular  cases,  steady-state  linearized  theory  also  in  general 
yields  infinite  pressures  at  3=0. 


rwy,T^r-^n>B  ■»-  - 


jyywjjS  jjfc  )j*;i  y-  iiHWWjWiWj 
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infinitely  long  fuselage  are  slightly  indented  in  the  vicinity  of  the 
wing  and  possess,  according  to  the  theory,  the  important  property  of 
having  zero  wave  drag. 


Velocity  potential  at  the  wing  surface.-  Let  the  surface 
zu(x,y)  =  -zi(xfy)  of  a  symmetrical  wing  be  specified  and  assume  that 
the  indented  fuselage  is  a  body  of  revolution  which  deviates  slightly 
from  a  basic  circular  cylinder  of  radius  Rq.  The  cross-sectional  area 
S(x)  of  the  exposed  wing  is  then  given  by 


J  Uo 


S(x) 


(x) 

Zu(x,Tl)dT| 


R(x) 


(29) 


where  R  is  the  radius  of  the  indented 
body.  (See  sketch  (e).)  In  order  for 
the  streamwise  gradient  of  cross- 
sectional  area  of  the  entire  combina¬ 
tion  to  be  zero,  the  relation 


«R2(x)  =  «R02  -  S(x),  S(x) <  <  «R02 

(30a) 


2«B  g  .  -SMx) 


(30b) 


must  hold.  Since  the  quantity  S(x)  is  known  from  formula  (29),  equa¬ 
tion  (30a)  can  thus  be  solved  for  the  radius  R  of  the  body. 


Now  the  perturbation  potential  in  the  plane  at  zi  =  0  has  the 


form 


1  fSl 

<Pa(x>yi>o)  mzl  vi(x»Tu)in  |yi  -  m  I  4m 

-Si 


(31a) 


or 


®i 

<P2(x,yi,o)  =  Wi(x,ru)in|yi  -  ml  am  +  kj  wi(x,m) inlyi2  -  m2l4m 
-n  ri  (31b) 


i  ' 
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with 


R2  od 

y±  -  y  +  ~y  >  rx=2R 


(32) 


The  term  corresponding  t6  the  g(x)  of  equations  (8),  which  is  to  be  added 
to  this  Vz,  is  zero  since  there  is  no  change  in  the  cross-sectional  area 
of  the  specified  shape  of  the  yz  plane. 


From  the  last  relation  in  equation  (19a)  it  follows  (upon  putting  R 
in  place  of  Ro)  that  in  the  interval  m2  <  ri2  tne  boundary  condition 
may  be  expressed 


vr  _  (dR/dx)U0  _  ^  R 

2  *ln  »  ’  2  *7i  -  (Wn)s 


(33) 


For  very  thin  wings ,  and  for  deviation  Rq  -  R  of  the  same  order  aB  the 
wing  thickness,  the  condition  in  the  interval  rx2  <  ilj2  <8^  is  approxi¬ 
mately 


w  Oq^z^/^x) 

1  -  r2/^2  1  -  R2/n2 


(3*0 


Making  use  of  equation  (33)>  one  can  then  write  the  first  term  on 
the  right  in  equation  (31b)  as 

rx 

^  J  w1(x,q1)ln  |yx  -  Ti1|dq1  =  U0R  ~  In  y,  y2  >  R2  (35) 

-*•1 

where,  from  equation  (30b),  the  quantity  R(dR/dx)  may  be  replaced  by 
-l/2n  times  the  gradient  of  cross  sectional  area  S*(x)  of  the  wing. 

In  the  second  integral  in  equation  (31b),  set 


=  n  + 


R2 


and  use  relations  ( 3I) .  The  final  form  for  the  velocity  potential  at  the 
upper  surface  of  the  wing  thus  becomes 

u0 

<P(x,y,0)  *  ~  2n  s'(x)ln  y  + 


£  f U*, 


q)ln 


(y2  -  Tia)(y*Tf  -  R4) 


y®n2 


dq, 


y  >  R 


p— 


(36) 
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Formula  (36)  for  the  velocity  potential  for  a  slender  combination 
whose  body  is  indented  according  to  the  area  rule  (ref.  12)  is  evidently 
independent  of  Mach  number  and  therefore  holds  (within  the  assumptions 
of  the  simplified  theory)  for  subsonic,  transonic,  and  supersonic  speeds. 
If  the  limits  of  integration  are  adjusted,  the  equation  may  also  be  used 
for  combinations  with  sweptback  wings. 

Pressure  coefficient  at  the  wing  surface.-  On  the  surface  of  the  wing 
the  pressure  coefficient  can  be  found  from  relation  (l4a),  that  is, 


where  the  potential  cp  is  obtained  from  equation  (36). 


APPLICATIONS 


The  formulas  which  were  given  the  preceeding  part  of  this  report  will 
now  be  applied  for  the  purpose  of  performing  detailed  calculations  of  the 
pressure  distribution  for  some  particular  nonlifting  combinations  having 
a  symmetrical  wing  with  constant  chord. 


Wing  Mounted  on  a  Body  of  Circular  Cross  Section 


Pressure  coefficient  on  the  surface  of  the  wing.-  The  wing  of  the 
combination4  considered  here  is  mounted  on  a  circular  cylinder  and  has 
an  upper  surface  defined  by 

zu(x,y)  =  T(c0m  -  s  +  y)v/(s  -  y)(c0m  -  s  +  y)  (38) 


with 


s  =  mx  +  R0 

zu(x,y)  =  zu(x,-y)  =  -zz(x,y)  -  -zz(x,-y) 


and 


T  8t0 
3m2C  0&/T 

4This  combination  will  be  referred  to  as  the  basic  combination. 


(39) 
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the  constant  t0  being  the  maximum  thickness.  The  profile  resembles  a 
Joukowski  section  and  is  the  same  for  all  values  of  y.  On  the  starboard 
side  of  the  wing,  the  slope  of  the  surface  is  given  by  the  relation 


=  Mx>y)  = 


v„(x,y) 


=  (c0m  -  4s  +  4y) 


s  +  y 


s  -  y 


From  this  equation  it  is  seen  that  the  slope  possesses  a  square-root 
infinity  at  the  leading  edge  and  is  zero  at  the  trailing  edge. 

The  gradient  of  cross-sectional  area  of  the  wing  is 
Si' (x)  =  Vnn3(c0  -  x)Vx(c0  -  x)  =  4m  zu(x,R0),  (0  <  x  <  cQ)  (^la) 


in  region  1  in  sketch  (f)  and  is 

S2'(x)  =0,  cQ  <  x  <  (s0  -  R0)/m 
in  region  2;  for  region  3, 

S3'(x)  =  -4m  ^(xjSq),  (sq  -  R0)/m  <  x  <  (bq  -  R0  +  cQm)/m 

Since  S"  is  not  continuous  across  j 

the  line  AB,  the  approximate  theory 
employed  here  will  give  unrealistic 
results  in  the  regions  of  the  tips. 

In  calculating  the  pressures  for  the 
wing-body  combination,  attention  will  f 
therefore  be  confined  to  regions  1  I 

and  2  only.  The  calculations  to  be  c,° 

given  in  subsonic  flow  will  be  based 
on  the  assumption  that  the  tips  are 
located  far  downstream  and  have  no 
effect  on  the  other  two  regions . 

UBe  of  the  equations  (40)  and 
(4l)  in  formulas  (28a)  and  (28b)  yields 
the  following  final  results  for  the 
supersonic  pressure  coefficient  on  the 
upper  surface  of  the  wing  for  regions 
1  and  2: 


(4lb) 


(41b) 


y~mx*R0-c0m 


Sketch  (f) 


Region  1: 


-2 T  r  ,  /  X  ,  .  v  ho  -  X  8xy(c0  -  x) 

3r)  =  — —  l6RDtan  1  / -  +  2m(c0  -  4x)  /  — z - In  - — j~z - —rr  + 

Pflup  v  cQ  -  x  v  x  PcoCy2  -  R02) 


F(x,y)  +  F(x,-y)  -  2F(x,0)  +  f(x,^  +  f(x,  -  ^)J, 


(0  <  x  <  cQ)  (42a) 


Region  2 : 

^sun =  ®2t  ffiRo +  -  co)  + 


ye  -  c0m 
- - -  + 


t 

/  s  -  ccm  +  y  /  .  ,  Rq2\  A 

(cQm  -  4b  -  4y )J - ZTJ +  VC°m  "  kB  +  4  y  jJ‘ 

(=.■  - 


b  -  cQm  -  Rpg/y 
a  -  R02/y 


Rp2/y 

,2/y 

(cQ  <  x,  B  -  c0m  <  y  <  x;  s  =  mx  +  R0) 

(42b) 


where 

F(x,rj)  =  (c0m  -  4b  -  4,^ 


[f- 


tan' 


-i  (cQ  -  x)(2b  +  i})  -  (nx  +  RqCo) 


2 Vx(c0  -  x)(s  +  tj) ( 8  -  cQm  +  n) 


Roco)  1 

m  +  n)  J 


(43a) 


when  cQm  -  a  -  rj  <  0,  and 
F(x,ti)  =  (cQm  -  4s  -  4,>yS= 


■  8_i_H 
n 


c0(n  +  R0) 

(c0  -  x)(2s  +  n)  -  (nx  +  R0c0)  +  >/x( c0  -  x) ( b  +  l)(c0m  -  s  -  tj) 

(43b) 


when  com  -  8  -  1  >  0 
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It  is  seen  from  equation  (26)  that  the  formulas  for  the  pressure 
coefficient  in  subsonic  flow  can  be  obtained  by  adding  two  terms  to  the 
above  equations  for  the  supersonic  pressure  coefficient.  Since  the  value 
of  S’  in  region  2  is  identically  zero,  the  terms  to  be  added  for  region 
1  become 

c0 

~  tti'Sxf  ^  =  Tm3t4x  '  3c0],  (x  <  c0)  (Wte) 

0 


and  for  region  2 

46  =  Tm3[^X  ‘  3c°  +  (c°  '  . 

o  L 

The  subsonic  pressure  coefficient  on  the  upper  surface  of  the  wing  can 
therefore  be  written  in  region  1  as 


(x  >  c0) 


(44b) 


CP„ub  =  cPeup  *  -  3co>  +  1“  5^]  ^ 

and  in  region  2  as 

c  =  c  +  vtPrlkx  -  3c0  +  (cQ  -  4 x)J  *  C°H  (45b) 

psub  psup  L  v  x  j 

where  c „  is  given  by  equations  (42). 

^s  up 

Formulas  (42)  and  (45),  which  give  the  pressure  coefficients  on  the 
wing  of  the  combination  indicated  in  sketch  (f),  will  be  plotted  and  the 
results  discussed  in  a  later  section.  An  obvious  result  noted  now  is 
that  in  region  2  the  formulas  do  not  depend  on  Mach  number  for  either 
supersonic  or  subsonic  flow,  but  that  the  subsonic  coefficient  in  this 
region  is  always  greater  than  the  supersonic  coefficient  because  the 
inequality 


m3r 


4x  -  3c0  + 


(c0 


(*6) 


holds  for  all  values  x  ^  c0 
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As  a  check,  let  us  now  consider  the  asymptotic  behavior  of  equa¬ 
tions  (38),  (42b),  and  (4*ib)  far  outboard  along  the  wing.  For  this  pur¬ 
pose,  we  introduce  a  change  in  coordinate  system  by  means  of  the  trans¬ 
formation 


x  =  xnsin  9  +  ynsin  6 
y  -  R0  =  -xncos  e  +  ynsin  0 


(47a) 


|  Uo 


(See  sketch  (g).)  Setting  tan  fl  =  m, 
one  can  then  write 

mxn  +  yn 


x  = 


y  -  Ro 


Vl  +  m2 
myn  -  *n 
«Vl  +  m2 


>(47b) 


mx  -  (y  -  Ro) 

’  ji  ♦ 

x  +  m(y-Ro) 

Yn  =  ./lTm5  J 


Mjj  =  m/^1  +  m2 

Un  =  U0m/Vl  +  5* 


Pn  =  1/n/1  +  m2 


(48) 


Equation  (38)  thus  becomes 


8t0  ,  n3/*/„  *1/* 


(cn  -  *,)  '  <*n> 
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which  agrees  with  the  approximate  equation  employed  in  two-dimensional 
section  theory  for  a  thin  Joukovski  base  profile  of  thickness  ratio 

]l  Vm 

-  .  (See  ref.  13.)  The  pressure  coefficient  for  such  a  two- 

3*/3cn 

dimensional  airfoil  is 


Cp  - 


Jffi.  =  - 


8t, 


^  Sxn  3V3cnPn 

or,  written  in  terms  of  the  xy  system, 


(49a) 


Cp  =  - 


8t 


3>/lc02 


£ —  [3c0m  -  4(mx  -  y  +  Rq)] 


(49b) 


This  equation  is  in  agreement  with  the  asymptotic  expression  of  the  pres¬ 
sure  coefficient  for  region  2  obtained  far  outboard  along  the  wing  from 
either  equation  (42b)  or  (45b). 

Pressure  coefficient  on  the  body.-  The  formula  for  determining  the 
pressure  coefficient  on  the  surface  of  the  body  is  (ref.  4) 


_2_  39  1  ( 1  /dR\2 

uQ  u^r*  vse/  1“  \dxy 


(50a) 


When  the  body  is  a  circular  cylinder,  this  relation  reduces  to 


(50b) 


It  is  apparent  that  application  of  equation  (50b)  requires  a  knowledge 
of  the  value  of  the  potential  9  in  space,  so  that  formulas  (28),  which 
hold  for  z  =  0,  must  be  ’’analytically  continued."  One  way  that  will  in 
effect  accomplish  this  is  as  follows: 

Set 


5  =  y  +  iz  =  Roe10 
#6*0)  = 


and  form  the  function 


<p(x,C,0)  +  i*r(x,£,0) 


(51) 
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where  <p(x,5,0)  is  taken  from  equations  (23)  and  (24),  with  the  function 
t  defined  by5 


s(x) 

t(x,$,0)  =  J  wu 

c 


(x,T|)dTJ 


(52) 


The  constructed  function  (51)  will  obviously  satisfy  Laplace’s  equation 
in  the  yz  plane;  and  it  can  be  verified  without  difficulty  that  the 
boundary  condition  vr  =  0  is  also  satisfied.  The  coefficient  on  the 
body  can  then  be  calculated  from  the  formula 


where  the  symbol  Re  means  that  the  real  part  is  to  be  taken. 


(53) 


For  region  1,  the  final  result  obtained  for  the  pressure  coefficient 
on  the  body  in  supersonic  flow  for  the  basic  combination  considered  in 
this  part  of  the  report  is 


16Rq  tan 


Cq  -  X 


2F(x,0)  +  F(x,Roe10)  + 


F(x,R0e‘ie)  +  F(x,-R0elw)  +  F(x,-R0e"iy)  + 


iB x 


Re 


2m(c0  -  4 x)J 


c0  -  x  8x(c0  -  x)e 
- = —  In - 


i© 


pco(e2i0  -  1)R0 


,  ia  I cQm  -  s  +  Roe10 

in(c0m  -  4s  +  4R0eiS)  / - rr - 

V  s  -  Roe16 


(*  £  cQ) 
(54a) 


where  the  function  F(x,q)  is  given  by  equations  (43).  The  result  found 
for  the  supersonic  pressure  coefficient  in  region  2  can  be  put  in  the  form 

5 An  arbitrary  real  function  E(x)  may,  in  general,  be  added  to  the 
right  side  of  equation  (52).  Such  a  function,  however,  would  in  no  way 
affect  the  real  part  of  d$/dx,  since  IE  is  purely  imaginary. 
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•sup 


m®T 


>y^ 


8R0  +  2m(4x  -  c0)  /---■  C°  +  2(4s  -  c0m)  J -j-2-  + 


(Com  .  4s  -  4R0ei9)  /S  +-R°ei9  ’  C°m 


s  +  R0e 


10 


(c0m  -  4s  -  4R0e"ie) 


ys  +  Roe“i0  -  ccm 

s  +  Roe“i0 


(c0m 


4b  +  4E0eie)  /S  ‘  R°ele  '  Co" 


s  -  R0e 


10 


+ 


(=o«  -  4b  ♦  4B0e-‘«)  A  -  R.°e'lfl  ' 

V  S  -  Rne“10 


U  >  Co)  (54b) 


Corresponding  formulae  for  the  pressure  coefficient  on  the  body  in  sub¬ 
sonic  flow  are  obtained  from  the  above  two  equations  by  adding  to  them 
the  terms  on  the  right  of  equations  (44a)  and  (44b),  respectively.  For 
any  fixed  value  of  y  in  region  2,  relations  (45b)  and  (46)  furnish  the 
inequality 


CPsub 


(55) 


which,  as  mentioned  before,  is  also  true  for  the  coefficients  on  the  wing. 


Wing  Mounted  on  an  Indented  Body 


The  constant-chord  wing  whose  thickness  distribution  on  the  starboard 
side  is  given  by 

zu(x,y)  =  T(c0m  -  s  +  y) V (s  -  y)(cQm  -  s  +  y),  R(x)  £  y  £  s 

s  =  mx  +  R0  >(56) 

8tq 


3m2co2  >/! 
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will  now  be  combined  with  an  indented  body  of  revolution  which  deviates 
slightly  from  a  basic  cylinder  in  the  manner  that  makes  the  local  cross- 


1  «■ 


sectional  area  of  the  configuration 
a  constant.  The  wing  has  the  same 
profile  along  EE  as  the  one  in  the 
previous  example.  (See  sketch  (h).) 
Since  the  gradient  of  cross-sectional 
area  of  the  combination  is  set  equal 
to  zero,  the  wave  drag,  according  to 
the  slender-body  theory,  is  also  zero. 

The  cross-sectional  area  of  the 
wing  in  region  2  is  equal  to  the  con¬ 
stant 

s 

s^x)  =  4  J  zu(x,n)dT)  =  Cg3^-—  (57) 

s-c0m 

and  the  radius  of  the  indented  body 
in  this  region  is  therefore  the  con¬ 
stant 


(58) 


Corresponding  equations  for  region  1  are 


rS  ,  .  T  f  o  r,  -1  c0®  “  2s  +  2Ri  , 

Sx(x)  =  4  J  zu(x,T|)dTi  =  —  |3®3c03c08  - + 

Ri 

2[2(s  -  Ri)(7c0m  -  4s  +  4rx)  -  Scq^W  (B  ”  Ri)(c0m  -  8  +  ^~i)}(59) 


and 


RiU)  =  Ro  I1  ~  »  O<R2<RX<R0 

v  rtRo 


(60) 


.'V  , -  : 
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Equation  (60),  however,  does  not  give  an  explicit  value  of  the  radius 
in  region  1  because  the  function  Sx  itself  involves  Rx.  Moreover, 
the  function  Sx  is  a  transcendental  function  of  Rx,  so  that  equa¬ 
tion  (60)  can  only  be  solved  by  graphical  or  other  methods  of  approxima¬ 
tion.  By  means  of  the  mean-value  theorem,  we  can  write  the  equation  in 
the  form 

Si(x)  =  S:i(x)  +  4zu(x,d)(R0  -  Ri)  (61) 

where  d  is  a  certain  value  in  the  interval  Rx  <  d  <  R0,  and 

2(l4xco  -  8X2  -  3co2)*y  x(c0  -  x)  +  3c0scos-1  —  J  (62) 

is  the  cross-sectional  area  of  the  wing  in  region  1  neglecting  the  addi¬ 
tional  area  exposed  by  the  indentation.  Since  the  deviation  R0  -  Rx  is 
assumed  to  be  of  the  same  order  as  the  wing  thickness,  neglect  of  the 
term  Uz^XjdHRo  -  Ri)  will  evidently  introduce  only  an  error  of  the 
second  order  in  thickness.  For  very  thin  wings,  the  radius  Rx  may  thus 
be  approximated  by 

Rl(x)  *  (63) 

which  is  now  an  explicit  function 
of  x. 


Sketch  (i)  indicates  the  varia¬ 
tion  of  the  radius  as  a  function  of 
x/c0  for  T0/mc0  equal  to  0.1  and 
values  of  R0/mc0  equal  to  0.5  and 

1.0. 

Relations  (l4b)  and  (53)  can  be 
immediately  employed  for  determining 
the  pressure  coefficient  on  the  sur¬ 
face  of  the  wing  and  on  the  body  in 
region  2,  because  the  portion  of  the  body  in  that  region  is  a  circular 
cylinder.  Use  of  these  relations  with  equation  (36),  (57)>  and  (58),  then 
finally  yields  for  the  coefficient  on  the  upper  surface  of  the  wing 


10 - - - 

JL 

CofTt 

°0  .5  iO  /5  20  2.5  3.0 

X/Co 

Sketch  (i) 


o  /  \  m3r 

sn(x)  =  *72” 


|2(4s  -  3c0m)  +  2(hs  -  c0m)  +  (cQ®  “4s  ’ 


Cp  =  mSTj^C^s  -  3c 


/  ,  R22n\  /s  -  cnm  -R22/y  (  .  ,  Rg2>\  /s  ~  c0m  +R22/y 

(c0m  -  us  +  4  \°m ’ 4s  ‘4  *  V  “V  -J 


(s  -  c0m  <  y  <  s)  (64) 


s 
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and  for  the  coefficient  on  the  body 


*p2=«2T 


2(4s  -  3c0m)  +  2(4e  -  c0m)  /  +  (c0m  -  4b  -  4r 2eid)J Lz£&m.+ 

V  V  8  -  R, 


8  -  R2e 


Mi!  + 

i0 


-iS»  /  s  -  cnm  +  RPe"i0  ,  ,  ,  „  1SX  /  s  -  c0m  -  R2e*0 

Qm  -  4a  -  4R2e  )  / - Q  ■„  +  ( c0m  -  4s  +  4R2e  )  / - 2 - ~— 

v  a  -  R2e-10  w  v  a  -  R2ei0 


(c( 


(65) 


These  equations  are  employed  for 


.  R0 

—  co  m  +  •/  m2 


Cn^r 


the  equality  sign  giving  the  value  of  x  at  the  trailing-edge  fuselage 
junction.  The  two  formulas  do  not  involve  Mach  number  and  apply  for  sub¬ 
sonic,  transonic,  and  supersonic  flight  speeds.  Comparison  of  formula 
(64)  with  equations  (42b)  and  (45b)  shows  that  for  a  fixed  y  in  region  2 
the  inequality 


CP  cPsub  ^  cPsup 


(66) 


is  satisfied  when  x  >  c0.  Values  calculated  in  region  2  for  the  subsonic 
and  supersonic  pressure  coefficients  along  a  section  on  the  surface  of  the 
wing  of  the  basic  combination  will  thus  be  less  than  the  values  obtained 
for  the  coefficient  3p  on  the  wing  surface  of  the  indented  combination. 

In  region  1,  the  body  is  not  cylindrical  but  it  is  found  that  the 
squared  terms  in  the  pressure  relation  (37)  may  be  neglected  since  they 
contribute  only  quantities  involving  the  second  and  higher  order  in  wing 

thickness  that  are  small  in  comparison  with  -  .  Even  the  first 

term  in  relation  (37)  gives  rise,  for  the  particular  combination  consid¬ 
ered,  to  some  small  terms  of  the  second  order.  Such  quantities,  however, 
are  also  found  to  be  negligible.  The  formulas  to  be  presented  here  for 
the  pressure  coefficients  in  region  1  will  therefore  contain  (like  those 
for  region  2)  only  terms  of  the  same  order  as  the  thickness  of  the  wing. 


4 
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NUMERICAL  RESULTS  AND  DISCUSSION 


Graphs  of  the  pressure  coefficients  for  a  wing  alone  and  for  the 
two  wing-body  combinations  in  sketches  (f)  and  (h)  are  shown  in  figures  1, 
2,  and  3*  Figures  1  and  2  give  plots  for  the  coefficients  along  several 
spanwise  stations  of  the  basic  combination  (R0/c0m  =  0.5)  for  subsonic 
and  supersonic  flow  at  mp  -  0.5,  while  figure  3  presents  plots  of  the 
coefficient  for  the  indented  wing-body  combination  which  are  independent 
of  Mach  number.  The  plot  in  figure  4  gives  the  variation  of  the  pressure 
on  the  body  in  the  plane  of  the  wing. 

Figures  1  and  2  show  that,  except  at  the  wing- fuselage  juncture,  the 
pressure  coefficients  have  a  finite  negative6  value  at  the  leading  edge 
and  increase  to  finite  positive  values  at  the  trailing  edge.  On  sections 
which  are  cut  by  the  plane  x  -  c0  (i.e.,  passing  through  the  trailing- 
edge  fuselage  juncture)  the  slopes  of  the  curves  are  discontinuous .  It 
is  also  apparent  that  the  effect  of  the  presence  of  the  body  on  the  coef¬ 
ficient  does  not  extend  very  far  downstream  beyond  the  trailing-edge 
fuselage  juncture;  the  difference  between  the  calculations  for  the  wing 
on  the  combination  and  those  for  the  wing  alone,  for  instance,  are  too 
small  to  show  up  in  the  plots  for  sections  more  than  one  chord  length 
from  the  body. 

In  figure  3,  the  pressure  coefficients  for  nearly  all  sections  on  the 
wing  of  the  indented  combination  are  very  close  to  the  curve  for  the  two- 
dimensional  wing.  The  discontinuity  in  the  slope  of  the  curve  that  was 
quite  noticeable  along  sections  such  as  BB  of  figures  1  and  2  is  far  less 
apparent  in  figure  3* 

The  graph  in  figure  5  is  the  function  which,  according  to  equa¬ 
tions  (27),  can  be  added  to  the  values  given  in  figures  1  and  2  for  the 
subsonic  or  supersonic  pressure  coefficient  along  the  section  AA  and  BB 
of  the  wing  to  yield  values  of  the  coefficient  at  other  Mach  numbers. 
(Sections  CC,  DD,  and  EE  lie  in  region  2  where  the  coefficients  do  not 

depend  on  Mach  number.)  For  example,  the  pressure  - - - — at  the 

(To/c  0)m 

leading  edge  along  BB  in  subsonic  flow  for  mg  =  0.4  is.  using  figures  1 
and  5  in  conjunction  with  equations  (27),  -4.1  -0.2  =  -4.3. 

Figure  6  shows  isobaric  charts  of  the  pressure  coefficients  for  the 
wing  without  body  in  subsonic  and  supersonic  flow  at  mp  =  0.5,  and  also 

6  The  occurrence  of  a  negative  pressure  at  the  nose  of  a  Joukowski- 
like  section  is  the  result  of  the  thin-airfoil  simplification.  A  more 
accurate  theory  would  show  a  small  region  of  positive  pressure  (i.e., 
a  stagnation  point  with  the  maximum  value  equal  to  the  impact  pressure 
of  the  component  of  stream  velocity  normal  to  the  edge) . 
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a  chart  giving  the  two-dimensional  results.  The  figures  6(a)  and  6(b) 
indicate  that  the  pattern  of  the  isobars  in  the  region  behind  the  trailing- 
edge  fuselage  juncture  is  essentially  the  same  as  the  two-dimensional  in 
figure  6(c),  but  that  in  the  region  near  and  upstream  of  the  juncture,  a 
marked  deviation  from  the  straight  isobars  of  the  two-dimensional  case  is 
evident . 

Isobaric  maps  are  also  shown  in  figure  7  for  the  basic  wing-body 
combination  in  subsonic  and  supersonic  flow,  and  for  the  indented  wing- 
body  combination.  Figures  7(a)  and  7(b)  illustrate  that  in  the  region 
downstream  from  the  trailing-edge  fuselage  juncture  (even  near  the  juncture 
itself)  the  isobars  are  not  much  different  from  those  for  the  wing  alone, 
and  that  the  body  therefore  has  little  effect  in  this  region.  In  the 
region  adjacent  to  the  body,  the  pattern  of  the  isobars  is  qualitatively 
similar  to  the  case  for  the  wing  alone  but  the  pressures  are  lower.  A 
remarkable  difference  between  the  chart  in  figure  7(c)  for  the  indented 
combination  and  those  in  figures  7(a)  and  7(b)  for  the  basic  combination 
is  that  the  isobaric  pattern  on  the  wing  for  the  indented  combination  is 
essentially  two-dimensional  over  practically  the  whole  wing. 

Examination  of  the  three  charts  in  figure  7  also  indicate  that  the 
maximum  negative  pressure  on  the  wing  occurs  at  the  leading  edge  near  the 
boundary  between  regions  1  and  2.  In  fact,  it  can  be  shown  that  the  max¬ 
imum  occurs  at  the  boundary.  In  view  of  the  inequality  (66),  indentation 
in  accordance  with  the  area  rule  reduces  the  maximum  perturbation  veloci¬ 
ties  on  the  wing. 


Ames  Aeronautical  Laboratory 
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references 


1.  Monk,  Max.  M.:  The  Aerodynamic  Forces  on  Airship  Hulls. 

NACA  Rep.  184,  1924. 

2.  Jones,  Robert  T.:  Properties  of  Low-Aspect-Ratio  Pointed  Wings  at 

Speeds  Below  and  Above  the  Speed  of  Sound.  NACA  Rep.  835,  1946. 
(Supersedes  NACA  TN  1032.) 

3.  Ward,  G.  N.:  Supersonic  Flow  Past  Slender  Pointed  Bodies.  Quart. 

Mech.  and  Appl.  Math.,  vol.  II,  pt.  1,  1949,  pp.  75"97» 

4.  Heaslet,  Max.  A.,  and  Lomax,  Harvard:  The  Calculation  of  Pressure 

on  Slender  Airplanes  in  Subsonic  and  Supersonic  Flow.  NACA  TN  2900, 
1953. 


32 


NACA  TN  3674 


5.  Heaslet,  Max.  A.,  and  Lomax,  Harvard:  Supersonic  and  Transonic 

Small  Perturbation  Theory.  Section  D.,  vol.  6,  of  High-Speed 
Aerodynamics  and  Jet  Propulsion.  Princeton  University  Press,  1954. 

6.  Adams,  Mac.  C.,  and  Sears,  William  R.:  Slender-Body  Theory  -  Review 

and  Extension.  Jour.  Aero.  Sci.,  vol.  20,  no.  2*,  Feb.,  1953 > 

PP .  85-98 . 

7.  Keune,  Friedrich:  Low  Aspect  Ratio  Wings  with  Small  Thickness  at 

Zero  Lift  in  Subsonic  and  Supersonic  Flow.  Kungl.  Tekniska 
Hogskolan,  Stockholm.  Institutionen  fur  Flygteknik.  Tech.  Note  21. 
June,  1952. 

8.  Oswatitsch,  K.,  and  Keune,  F.:  Nicht  angestellte  Korper  kleiner 

Spannweite  in  Unter-  und  Uberschallstromung.  Zeitschrift  fur 
Flugwissenschaften,  Bd.  1,  Heft  6,  Braunschweig,  Nov.  1953* 

9.  Harder,  Keith  C.  and  Klunker,  E.  B.:  On  Slender-Body  Theory  at 

Transonic  Speeds.  NACA  RM  L54A29a,  1954. 

10.  Lomax,  Harvard,  Heaslet,  Max.  A.,  and  Fuller,  Franklyn  B.:  Integrals 

and  Integral  Equations  in  Linearized  Wing  Theory.  NACA  Rep.  1054, 
1951.  (Supersedes  NACA  TN  2252.) 

11.  Lomax,  Harvard,  and  Byrd,  Paul  F.:  Theoretical  Aerodynamic  Charac¬ 

teristics  of  a  Family  of  Slender  Wing-Tail-Body  Combinations. 

NACA  TN  2554,  1951. 

12.  Whitcomb,  Richard  T. :  A  Study  of  the  Zero-Lift  Drag-Rise  Charac¬ 

teristics  of  Wing-Body  Combinations  Near  the  Speed  of  Sound. 

NACA  RM  L52H08,  1952. 

13.  Allen,  H.  Julian:  General  Theory  of  Airfoil  Sections  Having  Arbitrary 

Shape  or  Pressure  Distribution.  NACA  Rep.  833*  1945* 


Figure  1.-  Subsonic  pressure  coefficient  on  upper  surface  of  wing  of  the 

basic  combination;  m0  *  0.5. 
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Figure  5.-  Function  to  be  added  to  give  values  of  pressure  coefficient  for  other  Mach  numbers 


Figure  6.-  Isobar ic  Charts  for  the  wing  alone; 
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